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Abstract 

We evaluate the effective actions of super symmetric matrix models on fuzzy S"^ x S"^ up 
to the two loop level. Remarkably it turns out to be a consistent solution of HB matrix 
model. Based on the power counting and SUSY cancellation arguments, we can identify 
the 't Hooft coupling and large N scaling behavior of the effective actions to all orders. 
In the large N limit, the quantum corrections survive except in 2 dimensional limits. 
They are 0{N) and 0{Ni) for 4 and 6 dimensional spaces respectively. We argue that 
quantum effects single out 4 dimensionality among fuzzy homogeneous spaces. 
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1 Introduction 



Although string theory promises to tame quantum fluctuations of space-time, it also created 
deep questions of its own. One of such questions is to explain 4 dimensionality of space-time 
since the fundamental dimension in string theory is rather 10 or 11. Traditionally invisible 
dimensions are assumed to be compactified at Planck scale. More recently, it has been 
recognized that we could as well be living on 4 dimensional branes with large or even infinite 
extra dimensions. 

In any case, we certainly need to derive 4 dimensional gauge theory and gravitation from 
string theory. This problem can be compared to quark confinement problem in QCD. Any 
nonperturbative formulation of QCD must explain it. In this respect lattice gauge theory has 
been recognized as such since confinement is a natural phenomenon in the strong coupling 
limit. 

We believe that matrix models arc promising approach to investigate these nonpertur- 
bative questions in string theory [1][2]. Through them, string theory communicates with 
another promising idea of quantum space-time namely non- commutative geometry [3] [4]. In 
fact, non-commutative(NC) gauge theory is naturally obtained in matrix models with non- 
commutative backgrounds [5] [6]. The gauge invariant observables in NC gauge theory are 
the Wilson lines [7]. They play crucial roles to elucidate the gravitational aspects of NC 
gauge theory[8][9][10]. 

In the context of IIB matrix model, this question has been addressed from several different 
methods. We can cite branched polymer picture[ll], complex phase effects[12] and mean-field 
approximations [13] [14]. Although the results are encouraging thus far, it is fair to say that 
the problem is still far from settled. In particular the 4 dimensional gaussian distributions 
which have been obtained in the mean-field approximation are not realistic space-time yet. 

In this respect, we find fuzzy homogeneous spaces to be more attractive. We have success- 
fully constructed these spaces using IIB matrix models[18]. In the semiclassical limit, they 
reduce to Kahler manifolds up to dimension 6. Although we have not constructed experi- 
mentally favored de-Sitter space yet, they are closely related. Locally we obtain maximally 
supersymmetric Yang-Mills theory. Hence our model provides nonperturbative formulation 
of supersymmetric gauge theory as well if we can ignore gravitation [15] [16]. 

In this paper, we investigate IIB matrix model with (and without) Myers terms[17]. By 
introducing a Myers term, we can construct non-commutative gauge theory on fuzzy sphere 
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at classical level[19]. In our previous work, we have investigated quantum corrections of 
matrix models on fuzzy sphere up to the two loop level[20]. By modifying the Myers term, 
we can construct higher dimensional manifolds. In this paper, we investigate the simplest of 
such 4 dimensional manifolds, namely x S^. 

We compute the effective action up to the two loop level. Based on the power counting 
and SUSY cancellation arguments, we can identify the 't Hooft couphng and large N scahng 
behavior of the effective actions to all orders. In the large N hmit, the quantum corrections 
survive except in 2 dimensional hmits. They are generically 0{N) and O(iVt) for 4 and 6 
dimensional spaces respectively. 

These fuzzy homogeneous spaces are possible backgrounds in IIB matrix model as well. 
Although they are not classical solutions, they may minimize the effective action at quantum 
level. With these motivations, we also compute the effective action of IIB matrix model 
without Myers term around such backgrounds. We indeed find that fuzzy S"^ x 5*^ is a 
nontrivial solution of IIB matrix model at two loop level. In conjunction with our estimates 
of the large N scaling behavior of the quatum corrections, we argue that 4 dimensionality is 
singled out among fuzzy homogeneous spaces. 

The organization of this paper is as follows. In section 2, we investigate effective action of 
the deformed IIB matrix model whose classical solutions contain fuzzy S"^ x S"^. In section 3, 
we investigate the effective actions for fuzzy S"^ x S"^ in IIB matrix model itself. We conclude 
in section 4 with discussions. In Appendices A and B, we explain detailed calculations of 
two loop effective actions on fuzzy x S"^ with and without a Myers term. 

2 Effective actions in matrix models 

NC gauge theories on compact homogeneous spaces G/H can be constructed through matrix 
models. For this purpose, we may deform IIB matrix model as follows [18] 

SiiB ^ SiiB + '-U.pTr[A^, A,]A„ (2.1) 

where ff^i,p is the structure constant of a compact Lie group G. Since there are 10 Hermitian 
matrices in IIB matrix model, the number of the Lie generators of G cannot exceed 10 
in this construction. Within such a constraint, we can realize fuzzy Kahler manifolds up to 
dimension 6 such as = SU{2)/U{1),CP^ = SU (3) /U (2) [21] or CP^ = SO{5)/U{2)[22] 
as classical solutions of matrix models. 
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Since these models possess the translation invariance 



A^^A^ + c^, (2.2) 

and also 

V'^V' + e, (2.3) 

we remove these zero-modes by restricting ^4^ and ip to be traceless. 
The equation of motion is 

[A^, [A^, A,]] + if^p.[A^, Ap] = 0. (2.4) 

The nontrivial classical solutions are 

A'^^e, other AI^^O, (2.5) 

where i"'s satisfy the Lie algebra of G or its sub-group. We have investigated quantum 
corrections in supersymmetric matrix models on fuzzy 5*^ [20]. In this paper we extend our 
investigations to higher dimensional fuzzy homogeneous space. Although we investigate the 
simplest of such manifolds: S'^ x S'^, it may reveal generic features of matrix models on 
homogeneous spaces. 

In order to obtain NC gauge theory on fuzzy 5"^ x S^, we choose G — SU{2) x SU{2) 
with the following /^i/pi 

fixup = /w; (/^,i^,P) e (8,9,0), 
fi^up = fe^up] (/i,z/,p) e (1,2,3), 
other fl^^s = 0. (2.6) 

We investigate the following classical solutions: 

< = /Jm®!; (/^ = 8,9,0), 

< = /I® J/.; (/^ = 1,2,3), 

other A^J; = 0, (2.7) 

where and arc angular momentum operators. Wc further assume that and act 
on the n copies of spin li and I2 representations respectively with N = n(2/i + l)(2/2 + 1). 
In this paper, we always assume that k are large which implies n « N. These solutions 
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represent n coincident fuzzy S'^ x S*^. We refer to the constituent spaces as branes in this 
paper. In the hmit li = I2, our solution represents a 4 dimensional space while it reduces to 
a 2 dimensional space (S*^) in the I2 — limit. 

The reducibility implies that there are — 1 linearly independent Hermitian traceless 
matrices which commute with the classical solutions. They form the Lie algebra of SU{n). 
Its Cartan subalgebra represents the relative center of mass coordinates of the branes. They 
will be called as zero-modes in what follows. 

The classical action associated with this solution is 

-^-^N{h{h + l) + h{l2 + l)). (2.8) 



In the 4d limit (/i = /2), it becomes 



JlNh{h + l)r.-^N\ (2.9) 
In the 2d limit (^2 = 0), we obtain 

-^-^Nh{h + 1) ~ (2-10) 

Since it is minimized when ^2 = and n = 1, a single fuzzy S'^ with U{1) gauge group 
is classically favored. As for the reducible representations corresponding to the multiple 
branes, the classical action does not depend on their relative positions(zero-modes). 

In our calculation of the partition function, we divide out the following gauge volume of 
SU{N)/Z]sr by gauge fixing 



2^ TT— ^ ^ , ,, , 2.11 
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which appeared as the universal factor in [23]. We recall that the 'exact' free energy of IIB 
matrix model is as follows in this normalization [24] 

-ME A)- (2-12) 

n|Ar " 

Let us denote the bosonic and fermionic zero- modes as and With the presence of 

zero-modes, we integrate massive modes first to obtain the Wilsonian effective action which 
is a functional of zero-modes. At the one loop level, we obtain[ll]: 

-TrlogiP') - \TTlog ((P^ + ^F,.r''^)(^^)) , (2.13) 



where = A'^^ + X/^ and 

[e,X]=SX (2.14) 

We first estimate (2.13) in the coincident hmit where the bosonic zero-modes are smalL 
Since the leading contributions in the large limit come from large eigenvalues, we expand 
(2.13) into the power series of F^,^, and 5. The leading contribution is 



1 „ o^.„. 1 



= 2rr^ ^ 2„^ E(2, + l)(2p + l,-^--^^-^. (2.15) 

In the 4d limit (Zi = ^2), we evaluate it as 

Alog{2)nN. (2.16) 

In the 2d limit {I2 = 0), it becomes 

AnHog{N/n). (2.17) 

In this process, we also obtain the products of the following polynomials which contain 
fermionic zero-modes: 

I - 1 

S^,u = ^Sr^^TT^ri^S, (2.18) 

where m < S{n — 1). Since these terms are less singular than (2.15) in the large N limit, we 
may estimate TrS'^ ~ 0(1/ f^'^). After the fermion zero- mode integration, the normaliza- 
tion of the bosonic zero-mode integration measure is determines as 

In the presence of n coincident branes, we have found that (2.16) and (2.17) scales as 
with N/n being fixed. Since such configurations are further suppressed by the phase space 
of bosonic zero-modes, we conclude that the branes tend to move away from each other. 
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If two branes are separated by a distance which is much larger than their radii /, we 
can approximate TrS^i, as 

TrS,. - i^Sr,^r.S. (2.20) 

After integrating fermion zero-modes, we obtain the following potential between them with 
respect to the identical bosonic zero- mode integration measure as in (2.19). 

2Alog{\x\/l) + Slogl {2d limit), 

2Alog{\x\/l)-Slogl {Ad limit). (2.21) 



On the other hand, (2.15) can be estimated as 
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-m^ — . (2.22) 

x^ 

In the limit of \x\ » /, the former dominates the latter. Therefore the bosonic zero-mode 
integration converges at large distance and branes cannot move away from their neighbors 
much farther than their radii. 

Let us consider loosely bound branched polymer like configurations of branes which are 
separated from their neighbors by a distance d » I. The one loop level effective action 
corresponding to such a configuration can be estimated as 

mog{2)nf + (n - l)log{fyi^^) + U{n - l)log{d/l) {M limit), 

Anlog{2l) + {n-l)log{fH) + U{n-l)log{d/l) {2d limit). (2.23) 

This is the best upper bound of the one loop effective action we can obtain so far. We can 

interpret the first term in each limit as the one-loop self-energy of branes and the remaining 

terms as their interactions. We can trust the one loop estimation of brane interactions as 

long as they are well separated. Since the effective action is bounded by 0{n) quantity, it is 

consistent with our argument that n coincident branes cannot overlap each other. We thus 

argue that U{n) gauge symmetry is broken down to U{1) at the one loop level. 

We find it likely that branes settle into branched polymer like configurations by barely 

touching each other. ^ From the both limits in (2.23), we observe that the one loop effective 

action favors the 2d space {S^) over the 4d space (5*^ x S^) just like the tree action. 

^In the case of the 3d model with two component Majorana spinor [19], the fermion zero-mode integration 
results in the vanishing partition function. Such an effect may suppress the formation of the branched polymer 
like configurations. 
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We move on to study two loop corrections. Our strategy has been to obtain the Wilsonian 
effective action by integrating massive modes with fixed zero-modes. We delegate the detailed 
evaluation of the two loop effective actions to Appendix A. We first consider the 2d limit 
{I2 — 0) which has been studied in [20]. The only novelty here is that we have the twice 
contribution from the Myers term (diagram (c) in Figure 1). 

The two loop effective action with U (n) gauge group is 

F{1,0) = -l(n^(-40F|'(/) + 45F5(0)-n(-40F3"^(/) + 45F5(0)) 



where 



'^^ f^{2l + lf ^''' " + 0(1/Ar2). (2.24) 



2Z + 1' 

We summarize the effective action up to the two loop to the leading order oll/N os 

;N^ + 4n^log{N/n)-35-—y{n^-l). (2.26) 



24n2 »v / / y47v 

We still need to evaluate the higher order corrections beyond the two loop. Although 
we do not evaluate them explicitly in this paper, we can estimate their magnitude to all 
orders in perturbation theory based on our investigations thus far. In Feynman amplitudes 
of matrix models on S'^, the momentum integrations of field theory are replaced by finite 
series. Our important observation is that the large N limit of the series can be estimated by 
the power counting arguments in field theory. It is well known that there is no ultraviolet 
divergences beyond the two loop level in 2d gauge theory. We may then conclude that all 
series are convergent in the large N limit in matrix models beyond two loop. 

From this observation, we can estimate the i-th loop planar contribution to be 0{n^^^/ {fHy~^). 
It is because a single factor oi l/l arises at each order due the 6j symbols in the inter- 
action vertices. It thus appears that the i-th loop contribution is n^O((A^j;^)'^^) where 
= 47rn^/ if'^N) is the 't Hooft coupling of the commutative U (n) gauge theory. In this 
way we can estimate the planar part of the effective action to all orders as: 

71 



^^-nN^ +AnHog{N/n) +n^h2{\lM)- (2-27) 

DA 



where /i2(A^jy^) denotes a certain function of 



From these arguments, we find the effective action in the 2d hmit is 0{N'^) in the large 
hmit with A^J^^ being kept fixed where the tree action dominates the quantum corrections. 
However the coincident hmits may be exceptional configurations since we have argued that 
U{n) gauge group is broken down to U{1) by the dissociation process of branes at the one 
loop level. 

In the case of a single brane with U{1) gauge group, the two loop effective action is 
suppressed by another power of 1/N: 

1 948 
F(Z,0) = -^^j,{Fm-Fni))--YN-,- (2.28) 

It is because the two loop amplitude is convergent in our model and the theory becomes 
free (ordinary U{1) gauge theory with adjoint matter) in the infrared limit. In other words, 
the planar and nonplanar contributions cancel to the leading order oi 1/N. Since the gauge 
couphng of NC C/(l) gauge theory is = Sn/f'^N'^, it is natural to find O(A^) quantum 
corrections at two loop. 

We can summarize the effective action up to the two loop to the leading order oil/N as 

-Y^N^ + ^log{N)-^. (2.29) 

Since we expect the same pattern in higher orders, the structure of the effective action to 
all orders would be 

-^N + Alog(N) + h',(X'). (2.30) 

where /ig denotes another function of A^. 

In the absence of UV and IR contributions, we are left with supersymmetric U{1) NC 
gauge theory on in the large N limit. The absence of 0{N) quantum corrections in (2.30) 
is consistent with the vanishing quantum corrections in the flat space. We thus argue that 
the effective action for a single brane is 0{N) with A^ being fixed. In such a large N hmit, 
the theory is again classical as the tree action dominates quantum corrections. 

With the presence of n branes, we have argued that the branched polymer like con- 
figurations with U{1) gauge group are preferred at the one loop level. By assuming such 
configurations, the effective action for n branes can be bounded from the above as before: 

~ ;N + Anlog{N/n) + nh'^{\^) 



3A2" 

-7(n - l)log{\) - 6(n - l)log{l) + 14(n - l)log{d/l), (2.31) 
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where ~ r? j j'^N'^ . The first and second hnes in (2.31) correspond to the self-energy and 
interactions of branes respectively. 

Since the classical action dominates quantum corrections in (2.31), a single brane mini- 
mizes it with j^N"^ being fixed. Nevertheless the whole configurations with different numbers 
of branes energetically degenerate in the strong coupling limit. In such a situation, we argue 
that multi-brane configurations become dominant because of their large entropy. 

We have pointed out that NC gauge theory on a fuzzy S'^ may be related to 2d gravity [20] . 
It is well know that 2d supergravity with the central charge c > 1 is unstable against branched 
polymer formation. Since deformed IIB matrix model on S"^ with ?7(1) gauge group naively 
corresponds to c = 8, the formation of branched polymers in the strong coupling appears to 
be consistent with such a duality. 

We next investigate the 4d limit with /i = l^- Although the effective action could be 
quartically divergent by power counting, it is only quadratically divergent in this model. 
Since it is dominated by UV contributions, the effective action for U (n) gauge group is well 
approximated by the planar contributions as 



-8-{Fi{l,l) + 2F!{l,l)) -89.3-^, (2.32) 
where we have used the following numerically estimates: 

Fi{l,l) = 3.24 + 0(1/0, 

Fi{l,l) = 3.96 + 0(1//). (2.33) 

We can simply put n — 1 in (2.32) for U{1) gauge group in contrast to the 2d hmit. 
The total effective action up to the leading order oi 1/N is 

•' -N^ + 4log{2)nN-89.3—. (2.34) 
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We can estimate the i-th loop contribution to be 0(n*+^/^/(/^/^)*~^)) since we obtain a fac- 
tor of 1//^ from 6j symbols in the interaction vertices. Here we also assume that the leading 
contributions cancel due to SUSY. Under the assumption, the amplitude is quadratically 
divergent which results in an over all factor of . It thus appears that the i-th loop contri- 
bution is 0{nN{}^)^-^) where = {^-Kfr? is the 't Hooft couphng of U{n) gauge 
theory. In this way we can estimate the planar part of the effective action to all orders as: 

( - ^ + 4/0^(2) + h^{\^))nN. (2.35) 
9 



where /i4(A^) denotes a function of A^. Therefore we find that the 4d effective action is 0{N) 
in the large N hmit with being fixed. 

The effective action of loosely bound branched polymer like branes can be estimated as 
follows 

{-^^ + ^logi2) + hi\'))N 

-7{n - l)log{\) - 19(n - l)log{l) + 14(n - l)log{d/l), (2.36) 

where A^ = (47r)^n/(/'^A^). Here again the first and second line correspond to the self- 
energy and the interactions of branes respectively. We observe that the interactions are 
sub-dominant in comparison to the self-energy since the former is at most 0{nlog{l)). Thus 
the favored 4d configuration is such that it minimizes the self-energy with f^N being fixed. 
Since 't Hooft coupling A changes with n, we need to determine h^{}?) before answering such 
a question. We will estimate the strong coupling behavior of /i4(A^) in the next section. 

Nevertheless in the large scaling region with f^N fixed which is appropriate in 4d 
limit, the effective action is negative and scales as A^^ in 2d limit . Therefore in such a weak 
coupling regime, a single fuzzy S"^ always dominates in this model. The situation will be 
different in IIB matrix model which will be investigated in the next section. 

We can further determine the large N scaling behavior of the effective action in 6d case 
{S"^ X S'^ X S'^) in an analogous way. Such a space can be obtained as a classical solution 
of a matrix model by modifying the Myers term. The free energy with U (n) gauge group is 
estimated as 

-afW + bnY - + • • • , (2.37) 

where = n(2/ + 1)^. a, b, c are calculable numerical coefficients. The i-th loop contribution 
can be estimated by power counting arguments as 0{n^~^^l^/{f'^y~^l^~^). We have also 
assumed here that the the leading contributions cancel due to SUSY in this model. Therefore 
the loop expansion may make sense if we fix A^ = n/ fH. In this way we can estimate the 
effective action to all orders as: 

(^-± + b + h{X^))nH\ (2.38) 

where hQ{\^) denotes another function of A^. We conclude that the 6d effective action is 
0{N^) in the large N limit with f*N^ being fixed. 
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3 Effective action of IIB matrix model 



In this section, we investigate the effective action of IIB matrix model itself on fuzzy S'^ xS'^. 
Our investigation parallels with that in the previous section. The classical action of IIB 
matrix model for a configuration in (2.7) is 

~Tr[A^,A,][A^,A,] 

= ^7V(/i(/i + l) + /2(/2 + l)), (3.1) 

where N = n(2/i + l)(2/2 + !)• If we fix /, the 4 dimensional hmit (S*^ x S*^) with li = I2 is 
energetically favored over the two dimensional limit (S*^) with I2 = 0. We can also observe 
that larger gauge groups are favored. This preference is just opposite to the deformed IIB 
matrix model with the Myers term. However we still need to minimize the action (3.1) with 
respect to /. Since it is not stationary, we need to investigate higher order contributions to 
extract physical predictions. 

The leading term of the one loop effective action in the large N limit is 

-Tr({^')F^:.F,^F^,F,,) - 2Tr({^')F^,F^,F,,F^,) 

+lTr({^')F,,F,,F,,F,,) + ^Tr{{^')F,,F^,F^,F^,) 



-p2 /it'^ A"^^ -^^^ A^y I 4 V^p2 

3^2 (2j + l)(2p + 1) g^2 (2j + l)(2p + + if +p'{p+ if) 

|;(j(j + i)+p(p+i))2 ^ {j{j + i)+p{p + i)y 

i„2y- (2J + 1)(2P + 1) .32. 

i;{j{j + i)+p{p + iW' ^ ■ ' 

In 2d limit, it is estimated as 
In 4d limit, we can estimate it as 

-nHog{N/n). (3.4) 

The one loop corrections are found to be sub-leading since they are much smaller than those 
in the preceding section. 

Our remaining task is to determine the structure of the effective actions to all orders. In 
2d limit, we can show that the tree action dominates over quantum corrections by repeating 



11 



the same arguments in the previous section. It is because the Myers terms are soft in the 
sense that they do not alter power counting arguments. Since we can derive the identical 
long range interaction (2.21) as well, branes tend to form branched polymers. With fixed 
/, tree action favors larger n in contrast to the preceding section. Hence the multiple brane 
configurations are always favored. The difference is that / is now a parameter which describes 
the scale of field expectation values. We need to minimize the effective action with respect 
to / as well. Since the tree action is not stationary, we may conclude that 2d homogeneous 
spaces are after all not realized in IIB matrix model. 

In 4d case, we can repeat the same argument with the preceding section to show that 
the effective action is of the following form. 
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{^^ + h,{X'))nN, (3.5) 

where = (47r)^n^/(/^A^). h4{X^) is a function which describes the quantum corrections 
beyond one loop. Explicit evaluations of the two loop contribution in 4d limit are reported 
in Appendix B. Our calculation shows that 

h,{X') = 1-63^ + • • • . (3.6) 

Our next task is to minimize the effective action not only with respect to n but also with 
/ or equivalently A^. Due to the existence of nontrivial quantum corrections h4{\'^), A^ could 
be fixed as A^ by minimizing (3.5) independently of n. At the two loop level, we find 

where the effective action assumes the minimum: 2.55nN. Remarkably we find that 4d 
homogeneous spaces S'^ x S"^ with U (n) gauge group are solutions of IIB matrix model at 
two loop level as they minimize the effective action. Since the effective action is positive 
definite and O(n^), we further conclude that U{1) gauge group is favored in the coincident 
limit. The discovery of a consistent solution representing a realistic 4 dimensional space-time 
in IIB matrix model is the most important result of this paper. 

Although the existence of 4d homogeneous spaces in IIB matrix model has been shown 
only at the two loop level, we can argue that higher order corrections do not destabilize it. 
It is because the effective action grows in the both weak and strong coupling directions as 
we will argue shortly . Since the effective action is beheved to be bounded by (2.12), we 
may assume that 0{N) term is positive definite. With such an assumption, we can further 
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conclude that a single 4d homogeneous space (^^=1) minimizes the effective action in the 
coincident limit. 

In the presence of n identical 4d branes, we still need to investigate their zero-modes. 
We find it likely that branched polymer like configurations with U{1) gauge group will be 
favored again. It is because the effective action for such a configuration can be estimated as 
follows 

(S^ + ^,(P))7V, (3.8) 

where wc have neglected the interaction terms which are at most 0{nlog{l)). We thus find 
that the effective actions for different numbers of branes degenerate in 4d limit to the leading 
order of A^. This degeneracy could be lifted by the next leading 0(/ ~ y/N) terms in the 
effective action. We have thus presented a scenario in which a single 4d homogeneous space 
S'^ X S"^ with U{1) gauge group is reahzed in IIB matrix model at quantum level. Although 
we hope to understand the dynamics of the zero-modes in more detail, it is beyond the scope 
of this paper. 

In the remainder of this section, we argue that the homogeneous spaces are continuously 
connected to branched polymers. At the classical level, we have fuzzy homogeneous spaces 
made of quanta. The distance of the neighboring quanta is 0{fl^). If we fix the coupling 
of NC gauge theory g ~ 1/ f^l, the distance of the neighbors is 0(1/ g^). If we fix ~ ^/ Pl^ 
as in 6d or 2d with U{n) gauge group, neighboring points are separated by distances of 
0{l^/g^). They are homogeneously distributed on the respective manifolds. 

We can investigate the quantum effects on the distribution of eigenvalues as follows. Let 
us introduce the norm of A as 1^4 p = TrAA^. We can also split A — A'^ + a where A'^ and 
a denote the classical and quantum fields respectively. 

In a single sphere case, we can obtain the following estimates up to the one loop level: 

fj 

j^TrA^A, = fl%, + -j^^log{l)6^,, (3.9) 



where Sf^i, is projected in the sub-space in which resides. 
In 4d case, 

1 



^TrA^A, = fl%, + ^S^, = XN^S^, + ^S^,), (3.10) 
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where 5^^ assumes the value only in S"^ x S"^ sub-space. 
In 6d case, 

^TtA^A, = fr5^, + ^5^, = Xm{5^, + ^5^,), (3.11) 

where 5^1, assumes the value only in x S"^ x sub-space. 

A naive extrapolation of these weak coupling expressions suggests the following correla- 
tion functions in the strong coupling limit for the both 4d and 6d cases: 

^TrA^A, ~ XnU^,. (3.12) 

In the strong coupling region, we can no longer hope to relate A to the classical field A'^' as 
in the weak coupling region. A may be better defined as the parameter to describe the gauge 
invariant observables in (3.12). 

The gauge invariant correlator in (3.12) measures the distribution of quanta. AA^^ can be 
interpreted as their extension. This scaling behavior with respect to N reminds us branched 
polymers whose fractal dimension is 4 [11]. For branched polymers of N quanta which are 
separated by a distance d, the corresponding quantity is (dN)^. The scaling behavior in 
(3.12) is consistent with the branched polymers whose d ~ A^. 

The large A regime corresponds to the widely separated eigenvalue distribution of A. In 
such a regime, IIB matrix model is well approximated by the branched polymers and whose 
effective action can be estimated as 

UNlogiX^). (3.13) 

This is our prediction for the universal strong coupling behavior for the effective actions of 
IIB matrix model in homogeneous spaces. Our arguments here equally apply to the effective 
actions in the preceding section. 

In 2d, the quantum fluctuation in (3.9) appear to be smaller than the classical field even 
with large g ~ 1//^/ which is always smaller than / in the large limit. However we have 
argued in section 2 that branched polymers with large numbers of constituents will be formed 
in the strong coupling limit which is consistent with (3.12). 

In this section we have found that the effective action of IIB matrix model possesses the 
minimum with finite gauge coupling A in 4d among the simplest fuzzy homogeneous spaces. 



14 



The effective action favors 4 dimensionality since it is 0{N) in 4d wliile it is O^N^) in 6d. 
On tlie otlier liand, it does not possess 2 dimensional solutions like 5*^ since the theory is 
found to be classical in the large N limit. What is significant in our observation is that 
quantum fluctuations naturally select 4 dimensionality. Furthermore homogeneous spaces 
may be smoothly connected with the branched polymers in the strong coupling limit. 

4 Conclusions and Discussions 

In this paper, we have investigated two loop effective actions of supersymmetric matrix 
models on a 4 dimensional fuzzy manifold S"^ x S'^. We find it remarkable that fuzzy S"^ x S'^ 
turns out be a nontrivial solution of IIB matrix model at two loop level. Based on the 
power counting and SUSY cancellation arguments, we have identified the 't Hooft coupling 
and large scaling behavior of the effective actions to all orders. In the large limit, the 
quantum corrections survive except in 2 dimensional limits. They are generically 0{N) and 
0{N^) for 4 and 6 dimensional spaces respectively. These general arguments validate our 
solution beyond two loop level. 

Although we have investigated the simplest 4 dimensional manifold, we believe our results 
capture generic feature of matrix models on homogeneous spaces. It is because we have 
employed power counting and SUSY cancellation arguments which do not depend on detailed 
group structures of homogeneous spaces. Therefore we believe that 4 dimensionality is a 
generic prediction of IIB matrix model among fuzzy homogeneous spaces. 

Although 4 dimensional distributions are clearly favored in the mean field approximation, 
we believe we have made a substantial improvements from the optimization point of view 
as well. We have obtained O(A^) free energy while it is 0{N'^) in mean field approximation. 
However we still need to lower it down to the prediction in (2.12). We certainly hope to 
understand what kind of space-time achieves such a feat. 
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Appendix A 



In this appendix, we evaluate the two loop effective action of NC gauge theory on S'^ x S^. 
We consider NC gauge theory with U{1) gauge group in the context of a deformed IIB matrix 
model with a Myers term. It is straightforward to generalize our results to the case of U (n) 
gauge groups. 

The evaluation procedure parallels to that of NC gauge theory on S'^. There are 5 
diagrams to evaluate which are illustrated in Figure 1. (a),(b) and (c) represent contributions 
from gauge fields, (a) and (b) are of different topology while (c) involves the Myers type 
interaction, (d) involves ghost and (e) fermions respectively. 




(a) (b) (c) (d) (e) 

Figure 1: Feynman diagrams of 2 Loop corrections to the effective action 
We expand matrices in terms of the tensor product of matrix spherical harmonics: 

jmpq 

V' = X] i^jmpqYjm ® Ypq, (A.l) 
jmpq 

where 

P/. = Jm®1 (/x = 8,9,0), 
Pm = I^Jm = 1,2,3), 

other p'^s^Q. (A.2) 

The summations over j and p run up to j — 2li and p — 2I2 respectively where we assume 
N — (2Zi + 1) X (2Z2 + 1). We exclude the singlet state j = p = in the propagators in this 
paper. 
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We adopt the following representation of Yj^n- 

l-s f I j I 



(Yjm) 

where they are normalized as 



-ly 



-s m s 



2j + 1. 



(A.3) 



The cubic couplings of the matrix spherical harmonics can be evaluated as 



(A.4) 



X 



(-I)27(2ji + l)(2j2 + l)(2j3 + l) 

( ji h 33 \ / Ji J2 is 1 



\ mi 1712 TTls 

We refer to [25] for (3j) and {6j} symbols, 
bosonic propagators 



I I I 



(A.5) 



modes aj^^pg and ghost modes bjmpq, Cjmpq as follows 



From the quadratic terms in the gauge fixed action, we can read propagators of gauge boson 

as follows 

7 \ 7T Sjij2^piP2^mi-m2^qi-q2i 
-pl[Pl + i-) 

ni+qi 

(A.6) 



^jimipigi^j2m2P2q2 



( CjiTOipigi^j2m2P2g2 ) 

In terms of fields 



T' + + ^njA^P.5m,-m25q,-q2. 



b = 



c = 



^jmpq^jm <8 Ypg, 

impq 

^ ^ bjmpqYjm ® 
^ ] ^jmpqYjm ® 



(A.7) 



jmpq 



propagators become 



(cb) 



- E - 

1 (^—I'j'm+q 



(A.8) 
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contribution from 4-gauge boson vertex (a) 

The interaction vertex which involves 4 gauge bosons is 



/4 



It gives rises to the following contribution 



< "Yi^^ >lPI-21oop 

+ 1) + I))(j2(j2 + 1) +P2(P2 + 1)) 

Here we can use the completeness condition: 

We thus find 
45 1 

2 /■4 5-^ 5-^ 5Z 5Z 

iii2i3 mim2m3 P1P2P3 <?1'?2<?3 

(Tr Yj^jn^Yp^q^ [^jimi^piqi ■> ^j2m2^P2q2\) 

O'lO'l + 1) + 1))0202 + 1) +P2(P2 + 1)) 

- -45^ E E E E 

313233 ■>n\m2mz P1P2PZ ?i?2?3 

(2ji + l)(2j2 + l)(2j3 + l)(2pi + l)(2p2 + l)(2p3 + 1) 

(il(il + 1) + + I))(i2(i2 + 1) + P2(P2 + 1)) 

^ ( jl 32 js ^ / Pi P2 P3 W jl 32 js \ f ^2 P3 1 

\^ mi m2 ms J y qi q2 Q3 J \ h h h } \ k k k J 
= -45^ (Ff (Zi, k) - F^'ih, k)) , (A.12) 



(^l _ (^_iyi+32+33+Pl+P2+P3^ 



where 



(2ji + l)(2j2 + l)(2pi + l)(2p2 + l) 

+ 1) + + I))(j2(j2 + 1) +P2(P2 + 1)) ' 

jlj2 V\V2 

(2ii + l)(2j2+l)(2pi + l)(2p2 + l) f /l /l J2 1 f k k P2 



X 



(il(jl + l)+Pl(Pl + l))(i2(i2 + l)+P2(P2 + l)) Wl ^1 il j ^2 Pi }' 

(A.13) 



18 



Here we have used the property of 3j-symbol: 



y y ( j2 js Y = y = 1 (A.14) 



and the property of 6j-symboh 

E(2ji + l)(2p, + l)| 



jipi 



Jl J2 J3 [ i Pi P2 P3 Y 

1 h h h 1 



(2/i + l)(2i2 + l)' 

I k k js \\k k P3 ] ' 



Using the following relation, 



(A.15) 



contribution from 3-gauge boson vertices (b) 

Firstly we calculate 2-loop IPl contribution from 3-gauge boson vertices: 

V3 = yTr [p/,, a^][a^, a^] - yTr [p^, aj[a^, a^]. (A. 16) 

The result can be expressed in a compact form as: 

1 1/ P'^ - P ■ O P'^ — P ■ O \ 

{ ^VsVs ) = (10 - l)-^( « -^^Q^ »p - « p2g2^2 »np), (A.17) 

where P,Q,R are defined as 

Q/i^'2m2^2?2 = [P/i) ^■2m2^2?2]) 

-R/i^'ams^sgs = [P^i? ■^'a^is^sQs]- (A. 18) 

We have also introduced the following average: 

«X»p = y ^'*23^^'i23, 
ji,Pi,mi,qi 

«X»np = y ^IgsX^^isa, 

ji ,pi,nH,qi 

^123 = '^^^3li'ni^32m2^j3m3TrYp^q^p^q^Yp^q^. (A. 19) 



p2 _ p2 _ q2 

P ■ Q*123 = 7, ^*123, (A.20) 
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we evaluate (A. 17) as 



2/4 I « p2Q2 »P « p2Q2 » -P ) 

jljij-i mim2m3 piP2P3 yi'32y3 

(2ii + l)(2j2 + l)(2j3 + l)(2pi + l)(2p2 + l)(2p3 + 1) 



X 



(J2(j2 + 1) + P2{P2 + I))(j3(j3 + 1) + P3{P3 + 1)) 



^ ji 32 is y / pi P2 P3 y i ji 32 j3 y \ pi P2 P3 Y 

\mi 1712 rn^ ) \ qi ^3 ) \ k h h j \ k k k J 
^ (pfih, h) - F^ih, k) + Fiih, h) - F^^ih, k)) , (A.21) 



where 



Fiih, k) = F^ik, k) = T-, h ^ E Z^-^^^^^^^^^^,, . (A.22) 

Here we have used the following property of 6j-symbol: 



(2ji + l)(2pi + l) 

Jipi^OO 



\ 2 f \ 2 

jl J2 33 1 \ Pi P2 P3 1 



h h h \ 1 ^2 k l 



2 



1 '^i2j3'^p2,p3 



(2/i + 1) (2/2 + 1) (2/i + 1) (2/2 + 1) (2j2 + 1) (2p2 + 1) ' 

\h+puo^,^^\(Or..M^^^J^ J3\'fpi P2 P3 

k k k 



y: (-ir+^n2ji+i)(2pi+i) f f f 



^1 h 32 \ ] k k P2 \ ^j2,j3^p2,p3 



(A.23) 



[h h 33 j\k k P3 j (2/i + l)(2/2 + l)(2j2+l)(2p2+l)' 

It is because the singlet state is absent in the gluon propagator which results in the extra 
terms on the right-hand side of (A.23). Although they cancel each other in (A.21) for U{1) 
case, it is not the case for U (n) gauge groups. 

ghost contribution (d) 

We calculate the contribution from the ghost-gauge boson vertex: 

Vgn^-fTr\p„b][c,a^]. (A.24) 

The result is that: 
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contribution from cubic vertices (c) 

The cubic vertex which contains the structure constant of SU (2) x SU (2) is 

Vcubic = -^fUvpT^ a„]ap. (A.26) 



Their contribution is 



where 



< -;^VcubicVcubic >lPI-21oop 
8/1 1 



8 

7 



jFi'(/i,/2)-F3"^(/i,/2)), (A.27) 



(2ji + l)(2j2 + l)(2j3 + l)(2pi + l)(2p2 + l)(2p3 + 1) 

+ 1) + I))(j2(j2 + 1) +P2(P2 + I))(j3(j3 + 1) +P3(P3 + 1)) 

f Jl J2 J3 \ f Pi P2 P3 1 
1 /l /l h ] \ k k k ] ' 

F^^ih^k) = E E (-iyi+^'^+^'^+^i+p^+^^ 

ilj2,j3 Pl.P2,P3 

(2ii + l)(2i2 + l)(2i3 + l)(2pi + l)(2p2 + l)(2p3 + 1) 



X 



(jlO'l + 1) + l))(j202 + 1) +P2(P2 + 1))0303 + 1) +P3(P3 + 1)) 



fermionic contribution (e) 

The fermion propagator is 

11 1 1 



= J. (^r-^M + ^(^)V..pr^'^'^PpP. + 0{{^f)) . (A.29) 

We have expanded it in powers of 1/P^ up to the first nontrivial order which is sufficient to 
calculate the leading contribution to the two loop effective action. 
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2-loop (IPI) contribution to the effective action from the fermion-boson vertices is 

7^ V « P2g2^2 »P « (p2)2(g2)2^2 »P « (p2)2(g2)2/22 »f 

P-Q P^Q^ P-QP-Q 

^ PWP" ~ (P2)2(Q2)2p2 »"P « (p2)2(Q2)2p2 »«P 

= ^ [2Ff (ii, /2) - {h, h) + ^^(^1, k) - F^ih, k)] - (A.30) 



where 



E E 

JU'2J3 P1P2P3 

^ (2ji + l)(2j2 + l)(2j3 + l)(2pi + l)(2p2 + l)(2p3 + 1) 

+ 1) + I))2(j2(j2 + l)+P2{P2 + iWUsUs + 1) +P3(P3 + l)) 

X (j3(j3 + 1)P3(P3 + 1) + 2ji(ji + + 1) - 4j3(j3 + + 1) 

+2il(il + 1)P2(P2 + 1) + 4ji(ji + I)j2(i2 + 1)) 



X 



Jl J2 J3 I I Pi P2 P3 

h h ^1 I 1 h h h 



^_]^^ji+j2+j3+pi+P2+P3 

313233 P1P2P3 

(2ji + l)(2j2 + l)(2j3 + l)(2pi + l)(2p2 + l)(2p3 + 1) 



X 



+ l)+PliPl + I))'(j2(j2 + 1) +P2(P2 + I))'(j3(j3 + 1) +P3(P3 + 1)) 
X (j3(j3 + 1)P3(P3 + 1) + 2ji(ji + + 1) - 4j3(j3 + l)pi(pi + 1) 

+2jl(jl + 1)P2(P2 + 1) + 4ji(ji + l)j202 + 1)) 



2-loop effective action 

In this way we find the total 2-loop free energy h) of U{1) NC gauge theory on fuzzy 
5"^ X 5"^ as follows 

-F{h,k) = 8y^{Fi{h,k)-F7{h,k) + 2F!{h,k)-2Fr{h,l2))- (A.32) 
For U (n) gauge group, we obtain 

-F{h,l2) = j^(n\8Fi{h,h) + 16F!{h,l2)-^5F,{h,k)) 
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(A.33) 



We numerically find that planar contributions dominate over nonplanar contributions in 4d 
limit. 





= 3.24 + 0(1/0, 


F!{l,l) 


= 3.96 + 0(1//), 


Fs'iU) 


= 0(1//), 


Fr{l,l) 


= 0(1/0, 


Fm) 





(A.34) 

As the planar amplitudes are quadratically divergent, our findings are in accord with generic 
expectations in NC gauge theory. 



Appendix B 

We evaluate the two loop effective action of IIB matrix model in the 4d limit in this appendix. 
The principle difference caused by the introduction of the Myers term resides in the gauge 
boson propagators. The inverse gauge boson propagator in IIB matrix model with S"^ x S"^ 
background is 

P''5^r. + 2lf^,pP''. (B.l) 

Since the leading contribution to the two loop effective action comes from the quadratically 
divergent part, we can approximate the propagators as follows: 



5,.^^ - 2tU.,P''{^r + 4V(P)(-^)^ 
I,ji,v{P) are the following transverse tensor in S"^ x S'^: 

T (P) P"^ - P P + S P"^ - P P 



(B.2) 



(B.3) 



where and denote the components in the first and second three dimensional sub- 
spaces. The modification of the propagator is caused by a single or double insertions of 
—2iffMi^pP'^/P^ vertices into the minimal propagator. 

Strictly speaking there is a zero mode in the gauge boson propagator which corresponds 
to the total angular momentum 2 in the both SU{2). We may postpone to integrate this 
mode by considering the Wilsonian effective action. This problem does not modify the 
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leading 0{N) term of the effective action since it arises from the quadratically divergent UV 
contributions as we shall find. 

Since we do not have the Myers term in JIB matrix model, we should exclude the diagram 
(c) . We also need to consider the effects caused by the modification of gauge field propagators 
in (B.2). 

Prom diagram (a), we find the following extra planar contributions 

72 

-f,G, + y^G,, (B.4) 

where the first and the second hue in (B.4) is due to the insertions of —2if^^pPf/P'^ ver- 
tices into the identical and different propagators respectively. We also have introduced the 
following functions. 

Ci = ^ E (2ji + l)(2pi + l)(2j2 + l) (21^2 + 1) 

jl,Pl,j2,P2 

^^MJI + 1) + 1)^ J2(j2 + 1) +P2(P2 + 1)' 

G2 = «{^,n^,rR'»p- (B.5) 

From (b) and (d), we find 



-^(SG; + 6F3 + 32Hi - 32H2) 

+ 2G2 + 8F3 - 16Hs + I2H4 - 47/5), (B.6) 



where the first and the second hne in (B.6) is due to the insertions of —2ifn^pPf/P'^ vertices 
into the identical and different propagators respectively as before. ^ 
In this expression, we have introduced the following functions: 

H, = «P,PJAR){^f^2^,»P^ 



^In this evaluation, the following partial integration formula has been used. 



24 



H2 = «Pf^QulMi-B2f-5^7^»p^ 



1 1 



= «Q,Rj,.{p){—,y^{j^r»„ 



We can evaluate them as 

r> r m'^S- (2j + l)(2p+l) 



= E E (2ji + l)(2j2 + l)(2j3 + l)(2pi + l)(2p2 + l)(2p3 + l) 

(i3 03 + l)+P,3(P3 + l)) 



X 



+ 1) + 1))'(J2(J2 + 1) +P2(P2 + 1))^ 

h 32 h I f Pi P2 Pc 

l\ l\ W \ 1 I2 I2 I2 



Prom (e) 



-^(64G';-32F3 + 64i72) 

+-^(32G'; - I6G2 + UH^). (B.IO) 

Here the first hne in (B.IO) is due to the double insertions of —2if^vpPP/P'^ vertices into the 
gauge boson propagator and the second line is due to the modifications of the both gauge 
boson and fermion propagators. 

The total two loop effective action in IIB matrix model for 5"^ x S"^ space is 

-F = -^(16Gi- 8^2 -18F3 + 32i7i + 32/^2 + 48i/3 + 12i/4-4i^5 + 16F4) 

= -^(16G'i-8G'2-18F3 + 32//3 + 16F4), (B.ll) 

since Hi + H2 — 0, + H4 — and H3 — H^. The only remaining independent function 
we need to evaluate is H4: 

H ^ I \- \- (2ii + l)(2j2 + l)(2j3 + l)(2pi + l)(2p2 + l)(2p3 + 1) 

2 ,,£3p,i:p3 + 1) +PM + WihU2 + 1) +P2(P2 + I))^(i3(j3 + 1) +P3(P3 + 1)) 

X ( - 2(ji(ii + 1))^ - 03^3 + 1))' + 2il(jl + l)i202 + 1) + ^jlUl + + 1)) 



2 



^We can further show that F = 2^3//^. 
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We have numerically estimated the following functions in the large / limit as 

2Gi(/, /) - G2{1, 1) = 1.93 + 0(1//), 
H^{Ll) = 0.822 + 0(1//), 

F(/,/) = -^6.53 + 0(1//). (B.13) 
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